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Abstract. The impact of a kination-dominated (KD) phase generated by a quintessential expo- 
nential model on the thermal abundance of Supersymmetric (SUSY) extremely Weekly Interacting 
Massive Particles (e-WIMPs) is investigated. For values of the quintessential energy-density pa- 
rameter at the eve of nucleosynthesis close to its upper bound, we find that (i) the gravitino (G) 
constraint is totally evaded for unstable G's; (ii) the thermal abundance of stable G is not sufficient 
to account for the cold dark matte r (CDM) of the universe; (iii) the thermal abundance of axinos (a) 
can satisfy the CDM constraint for values of the initial ("reheating") temperature well above those 
required in the standard cosmology (SC). 
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INTRODUCTION 

A plethora of data [1] indicates that the two major components of the universe are CDM 
and Dark Energy (DE). The DE component can be explained with the introduction of 
a slowly evolving scalar field, q, called quintessence [2]. An open possibility in this 
scenario is the existence of an early KD era [3], where the universe is dominated by 
the kinetic energy of q. During this era, the expansion rate of the universe increases with 
respect to (w.r.t.) its value in SC. As a consequence, the relic abundance of WIMPs (e.g., 
the lightest neutralino) is also significantly enhanced [4, 5]. 

WIMPs are the most natural candidates for the second major component of the 
universe, CDM. In addition, supersymmetric theories predict the existence of even more 
weakly interacting massive particles, known as e-WIMPs [6]. These are the gravitino and 
the axino (SUSY partners of the graviton and the axion respectively). Their interaction 
rates are suppressed by the reduced Planck scale, mp = Mp/v^8^ (where Mp = 1.22 x 
10^^ GeV is the Planck mass) in the case of G and by the axion decay constant, 
fa ~ (10^^ — 10^^) GeV in the case of d. Due to the weakness of their interactions, e- 
WIMPs depart from chemical equilibrium very early and their relic density is diluted by 
primordial inflation. However, they can be reproduced in two ways: (i) in the thermal 
bath, through scatterings and decays involving superpartners [7, 8, 9], and (ii) non- 
thermally from the decay of the next-to-lightest supersymmetric particle; in this case 
the results are highly model dependent. 

In this talk, which is based on Ref. [10], we reconsider the creation of a quintessential 
kination scenario (QKS) in the context of the exponential quintessential model [5, 11], 
and investigate its impact on the thermal production of e-WIMPs. 



THE QUINTESSENTIAL EXPONENTIAL MODEL 



The quintessence field, q, of our model satisfies the equation: 

q + 3Hq + dV/dq = 0, where V = Vq^"^'^/'"? (1) 

is the adopted potential, dot denotes derivative w.r.t the cosmic time t and H is the 
Hubble parameter, H ~ y/pq + p^/ V^m^ where = q^/2 + V and = n^g.^.T'^ /30 the 
radiation energy density with g^^ the effective number of massless degrees of freedom. 
We impose on our quintessential model the following constraints: 

1. Initial Domination ofKination. We focus our attention in the range of parameters 
with Q}q = Q.q{Ti) > 0.5 where Q.q Pq/{pq + p^C) is the quintessential energy- 
density parameter. 

2. Nucleosynthesis (NS) Constraint. At the onset of NS, T^s = 1 MeV, pq is to be 
sufficiently suppressed w.r.t p^, i.e., [12] = Q^(^^fs) < 0.21 at 95% c.l. 

3. Inflationary Constraint. Assuming that the power spectrum of the curvature per- 
turbations is generated by an early inflationary scale, an upper bound on the initial 
value of H, Hi, can be obtained, namely Hi < 2.65 x 10^"^ GeV. 

4. Cosmic Coincidence Constraint. The present value of Pq, Pq, must be compatible 
with the preferred range for DE, implying Q.^^ = 0.74. 

5. Acceleration Constraint. Successful quintessence has to account for the present- 
day acceleration of the universe, i.e. —1 < w^(0) < —0.86 (95% c.l.), where 
Wq = (q^/2 — V)/ {q^/2 + V) is the barotropic index of the ^-field. 

Solving Eq. (1) with q{Ti) = and q{Ti) such that the condition 1 is satisfied, we find 
that, during its evolution, q undergoes three phases (see Fig. l-(a), where we plot logp,- 
with / = ^ and R versus T for Ti = lO'^ GeV, Q^^ = 0.01 and X = 0.5): 

• The KD phase [3], where Pq ~ q'^/2, implying ~ 1 and thus, pq T^. The 
transition from kination to radiation occurs at Tkr such that pq{T^R) = PR(7kR) - 

see Fig. l-(a). Solving this we find Tkr = T^s {g^Vgf^Y^^ ((1 - ''^ ■ 

• The frozen fleld dominated phase, where the universe becomes radiation dominated 
and pq is dominated initially by 4/2 and then by V. 

• The late-time attractor dominated phase during which p^ ~ V dominates the uni- 
versal evolution, with Wq ~ \Vq = A^/3 — 1 for A < a/3. 

Today we obtain a transition from the frozen-field to the attractor dominated phase. 
Although this does not provide a satisfactory resolution of the coincidence problem, the 
observational data can be reproduced. In particular, satisfying the condition 5 implies 
A < 0.9 whereas the condition 4 can be fulfilled by conveniently adjusting Vo [5, 11]. 
These two constraints are independent of the parameters Ti and which can be 
restricted by the requirements 1 - 3. In Fig. l-(b), we present the allowed region of 
our model in the log Ti — logH^^ plane (shaded in gray and light gray). We observe that 
for a reasonable set of the parameters (A,7i,n^^), our model can become consistent 
[5, 11] with the observational data. 
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FIGURE 1. (a) The evolution of logp; (p,- = pi/p^ with p^ = 8.099 x IQ-^^/jZ QeV^ andh = 0.72) 
with i = q (black line) and R (Ughtgray Une) as a function ofT for A = 0.5, Tj = 10^ GeV and Q^^ = 0.01 . 
(b) Allowed (gray and lightly gray shaded) region by the conditions 1 - 3 in the log 7i — log plane. 



THERMAL PRODUCTION OF SUSY e-WIMPs 

The number density of a SUSY e-WIMP x (where x stands for G or d) satisfies the 
Boltzmann equation [8, 9, 13] which can be written as 

nx + 3Hnx = Cxn''^^ + Y^:^mjTKi{mi/T)ri. (2) 

Here n^^ = ^{3)T^/n^ is the equilibrium number density of the bosonic relativistic 
species, m,- [g,] is the mass [number of degrees of freedom] of the particle / and 
is the modified Bessel function of the 2nd kind. In the relativistic regime (T ^ mi) 
has been calculated [7, 9] using the Hard Thermal Loop Approximation, resulting to 
= Cj'^, where 

^HT^ f (3;r/16C(3)mi)l3^i {l+Ml/3ml) CaglHka/ g a) for x = G, 
l(274/;rV^C(3))ln(1.108/g3) for X = a. 

Here, ga and Ma (with a = 1,2, 3) are the gauge coupling constants and gaugino 
masses respectively, associated with the gauge groups U{1)y, SU{2)i^ and 5f/(3)c, 
(ka) = (1.634,1.312,1.271) and (ca) = (33/5,27,72). Throughout our analysis we 
impose universal initial conditions for the gaugino masses, MaiMQUj) = M1/2 and 
gauge coupling constant unification, i.e., gaiMovr) — gGVT- Eq. (3) gives meaningful 
results only for T > Tc= 10^ GeV. Towards lower values of T, non-relativistic (T <S 
nii) contributions start playing an important role. In the case of d, where Q-ah^ takes 
cosmologically interesting values for T <^ nij, Q has been calculated numerically in 
Ref. [10]. In the latter case, F,- with i — g, q and B are taken into account, also, using the 
formulas ofRefs. [13, 14]. 

The impact of the KD phase on the relic density ^.^h^ = 2.748 x 10^ 7° m^/Ge\f 
can be found by solving Eqs. (1) and (2) numerically [10]. However, we can get a clear 




10'' 10"" 10' 1 10 itf lO'' 10"" 10' 1 10 itf 



mg(GeV) m^(GeV) 

(a) (b) 

FIGURE 2. Si-x^^ ^ function of (X = G [x = a]) for various Q.^^ 's, indicated on the curves, 
Ti = 10^ GeV and Ml/2 = 0.7 TeV [fa = lO" GeVJ (a [h]). ForQ^^ > 10'^^ we take in our computation 
the values ofrm indicated in Eq. (5). The CDM bounds are also, depicted by the two thin lines. 



picture of the results via semi-analytical estimates [10]. In the high T regime (Ti ^ Tq 
and Tkr » 7c) we find, within a 10% accuracy, that: 

yO-^r^fyo^'Ti in these, 

where s is the entropy density and)^^^ = /45 TtmpY^'^^C^^ with Y^'^ = rf^/s. From 

these expressions we can easily deduce that in QKS, decreases with respect to its 
value in the SC (being proportional to Ti in SC and to Tkr - lower than Ti - in QKS). 

In the low T regime {Ti <g Tq or Tkr <C 7c), we find cosmologically interesting 
solutions only in the case of a. Focusing on the most intriguing possibility, Ti ^ TsusY = 
1 TeV but Tkr <^ 7susY, and for the benchmark values of m,- used in our analysis: 

mq=\ TeV, mg= 1.5 TeV and m§ = 0.3 TeV. (5) 

we can write simple empirical relations which reproduce rather accurately the numeri- 
cal results. In particular, in the SC, using fitting technics, we get a relation with a 15% 

accuracy, namely: ^ah^ = Am« (1 -f-CTi) e'^/'^^/fl with A = 1.44 x 10^4 GeV, B = 
1 AS. All GeV and C = 0.001 /GeV. We observe that ^ah^ decreases sharply as Ti de- 
creases due to the exponential factor. In the QKS this suppression is avoided since 

Q.ah^ =Dma/fa \/^q^ ^ith D = 9.26 x 10^^ GeV. This relation reproduces the numer- 
ical results with excellent accuracy. We observe that Q.ah^ °^\/ ■\J^^^ or Q.ah^ <x Tkr. 

In Figs. 2-(a) [2-(b)], we display Cigh^ [Q.ah^] versus Mq [nia] for M1/2 = 0.7 TeV 
[fa — 10^^ GeV], Ti = 10^ GeV and various fl^^'s indicated in the curves. We observe 
that ^x^^ decreases as increases, D.Qh^ oc l/m^ and D.ah^ oc mg,. 
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FIGURE 3. In the logTi — logil^^ plane we show the (gray and lightly gray shaded) area, allowed 
by the quintessential requirements 1-3, and (a) the (black lined) area allowed by the G-constraint for 
= Ml/2 = 0.5 TeV; (b) the (black lined) area, allowed by Eq. (6), for G-CDM with = 0.1 TeV 
and 0.5 < Mi/2/TeV < 1 or (white lined) area, allowed byEq. (6), ford-CDM with ma — 5 GeV, w, 's of 
Eq. (5) and IQi" < fa/GeV < 10^^. 



QUINTESSENTIAL KINATION AND G-CONSTRAINT 

Unstable G can decay after the onset of NS, affecting by unaccepted amounts the 
primordial abundances of light elements. To avoid this, an upper bound on y^(rNs) 
can be extracted as a function of Mq assuming that the hadronic branching ratio of G 
is tiny. E.g., for iuq ~ 0.6 TeV, we obtain [15] Yg{TMs) ^ 10^^^ which means that 
^1^4x10^ GeV in the SC. This very restrictive upper bound on Ti can be avoided in 
the QKS, where we can set Ti= 10^ GeV. The upper bound on yg(rNs) can be satisfied 
for Q^^ > 10-21 or Tkr < 6.8 x 10^ GeV. The importance of a KD era in evading the G- 
constraint can be induced, also, by Fig. 3-(a), where we show the allowed (black lined) 
area by the G-constraint in the logTi — log^l^^ plane, for iuq — M1/2 — 500 GeV. The 
(gray and lightly gray shaded) area allowed by the quintessential requirements 1 - 3 is 
also shown. We clearly see that the G-constraint can be totally eluded in the QKS even 
with tiny values of almost independently on Ti. 

QUINTESSENTIAL KINATION AND G OR a CDM 

Stable x'^ constitute good CDM candidates provided that their relic density ^^^^ 
satisfies the CDM constraint [1] 

0.097 <%/j2< 0.12. (6) 

This constraint can be satisfied by both the G and d thermal abundance. However, in the 
case of G, Q.^^ is to be tuned to extremely low values whereas in the case of d, 
may be even close to its upper bound posed by the condition 2. Indeed, in Fig. 3-(b), 



we present the region in tlie logT — log^l^ plane allowed by both the quintessential 
requirements, 1-3, (gray and lightly gray shaded area) and Eq. (6) for G-CDM (black 
lined region) with thq = 100 GeV and 0.5 < Mi/j/TeV < 1 or a-CDM (white lined 
region) with ma = 5 GeV and 10^^ < /a/GeV < lO^^. Obviously Eq. (6) is met for a- 
CDM with much more natural ^2^^'s than those required for G-CDM. Therefore, d is 
more natural CDM candidate than G in the QKS. 



CONCLUSIONS 

We examined the impact of a KD epoch, generated by an quintessential exponential 
model, to the thermal abundance of G and d. The parameters of the quintessential model 
(A, Ti, O^^) were confined so as 0.5 < 0.q{Ti) < 1 and were constrained by using current 
observational data originating from NS, the acceleration of the universe, the inflationary 
scale and the DE density parameter. We found that < A < 0.9 and studied the allowed 
region in the (Tj, ^2^^)-plane. For unstable G, the G-constraint poses a lower bound on 

(almost independent of Ti). The CDM constraint can be satisfied by the G thermal 
abundance for extremely low H^^'s. On the contrary, this constraint can be fulfilled by 
the d thermal abundance with much larger H^^'s, making d a very good CDM candidate. 
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